Exact ground states for two new spin-1 quantum chains, 
new features of matrix product states 



S. Alipoui0, V. Karimipour H, L. Memarzadetjfl, 



Department of Physics, Sharif University of Technology, 
P.O. Box 11155-9161, Tehran, Iran 



Abstract 

We use the matrix product formalism to find exact ground states of two new 
spin-1 quantum chains with nearest neighbor interactions. One of the models, 
model I, describes a one-parameter family of quantum chains for which the ground 
state can be found exactly. In certain limit of the parameter, the Hamiltonian 
turns into the interesting case H = ^2,A>$i ■ S^+i) 2 . The other model which we 
label as model II, corresponds to a family of solvable three-state vertex models 
on square lattices. The ground state of this model is highly degenerate and the 
matrix product states is a generating state of such degenerate states. The simple 
structure of the matrix product state allows us to determine the properties of 
degenerate states which are otherwise difficult to determine. For both models we 
find exact expressions for correlation functions. 
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1 Introduction 



Quantum information theory and condensed matter physics, study many body sys- 
tems on lattices from complementary points of view. While in condensed matter 
physics, one starts from a Hamiltonian and seeks to determine the ground state 
in some approximate form, in quantum information theory the emphasis is on the 
properties of quantum states, for the quantification of which many tools have been 
developed in recent years. The subject of Matrix Product States (MPS) lies at the 
borderline of these two disciplines, since in this formalism, one starts from a quantum 
many body state with pre-determined properties, and then constructs a family of 
local Hamiltonians for which this state is an exact ground state. In this way, one may 
find interesting many-body systems for which the ground state and all its properties, 
i.e. correlation functions, can be calculated in closed analytical form, a very rare 
situation which is usually welcomed in condensed matter and statistical physics. 

The subject of MPS has a long history in condensed matter physics, the origins 
of which can be traced back to the work of Majumdar-Ghosh models p] which in 
turn inspired the construction of a larger family of solvable spin systems by Affleck, 
Kennedy, Lieb and Tasaki ( AKLT) in pj . The AKLT construction was further devel- 
oped in [3l H] under the name of finitely correlated spin chains or in [5j [6] under the 
name of optimal ground states. In its simplest version, which applies to translational- 
invariant systems on rings of N— sites, a matrix product state generalizes a product 
uncorrelated state by replacing numbers by matrices in the following way 

d 

= ^ tr(A il A i2 ---Ai N )\i 1 ,i 2 ,---iN}, (1) 

where Ai, i = 1, • • • d are a set of matrices, assigned to the local states |1), • • • \d) of a 
site. The normalization of these states is given by 

<V#> = tr(E N ), (2) 

where E := Ylt=i A* <8> A{. The dimensions of these matrices are arbitrary and are 
constrained by symmetry considerations and the details of model construction, i.e. 
the range of interaction. One can collect all the matrices in a vector- valued matrix A 
as follows 

d 

A = '£A i \i}, (3) 

i=i 

and write the matrix product state ([!]) as 

\ip)=tr(A®A---®A), (4) 

where the trace is taken over the matrix indices and the tensor product acts on basis 
vectors, i.e. 

tr(A <g> B) := Y,( A i)aA B j )p,<x\h j)- 
hi 
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The simple structure of the state @ allows an exact calculation of correlation 
functions. For example one and two-point functions of local operators are given by 



\0\ij)) tr(E E 



(V#) tr{E N ) ' 

MgfcgtM tr{E k - 1 E E l - k - l E E N - 1 ) 

(V#> tr{E N ) ' 1 ' 

where 

E = Y. A *i®A 3 (i\0\j). (6) 

In the thermodynamic limit N — > oo, the right hand sides of the above equations sim- 
plify even further, since in this limit, the eigenvalue(s) of E with largest magnitude 
dominate the traces. 

In recent years, this formalism has been used in developing exactly solvable mod- 
els in spin chains El El [TO1 EE] , spin ladders p21 031 O 02] , spin systems 
on two dimensional lattices [HI CE7J [18] , and the study of entanglement properties of 
spin systems near the points of quantum phase transitions [191 [20[ [21"| [22] . It has 
also been used extensively to find the stationary states of many types of stochastic 
systems of interacting particles in one dimensional chains, see for example [23, 24, 25J. 

A basic question is then whether we can construct general MPS and its parent 
Hamiltonians having a set of specific symmetries for quantum chains of spins. For 
spin-one systems, the first model was given by Affleck, Kennedy, Lieb, and Tasaki 
in [2], (not within the MPS formalism) which had full rotational symmetry, and was 
shown later to correspond to the matrix 

_ / |0) -V2|l) \ 

•A-aklt - /oiT\ ln\ > V'J 



-|0) 

where |1) = | — 1) with the parent Hamiltonian 

H = Sj ■ Sj+i + -(Sj • Sj + i) 2 . (8) 

i 

Then it was shown [5] that if one demands only rotational symmetry around the z 
axis in spin space, in addition to parity and spin-flip symmetries, a more general 
model can be constructed which is described by the matrix 



where g is a continuous parameter and a = ±1 is a discrete parameter. 

At first sight, construction of a matrix product state, finding its parent Hamil- 
tonian, and calculating the correlation functions, seems a straightforward procedure. 
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However if one demands symmetry properties, and more importantly demands that 
the final Hamiltonian have a physically interesting interpretation, then the problem 
will be quite non-trivial and interesting. Specially if one puts the formal MPS un- 
der scrutiny, one may be able to find many more states which are not themselves 
MPS representable, but have been captured by a single MPS in a nice way, i.e. as 
their generating state. This is what we will find for the models constructed in this 
paper. We believe that the richness of matrix product formalism has yet to be un- 
raveled by studying more and more examples. In this paper we try to construct two 
other spin-1 matrix product models, which have not been reported in the literature 
of matrix product states. The first model is a one-parameter family which has the 
interesting property to interpolate between two limits, namely between the Ising-like 
Hamiltonian 

tf^E^-iX^i+i- 1 ); ( 10 ) 

i 

and 

H 2 = £(S< • S m ) 2 . (11) 

i 

The ground state of (jlip . as we will see, breaks the rotational symmetry of the 
Hamiltonian. This is an example of the richness of matrix product formalism, that 
is by searching the space of solutions, one may come to corners where there are very 
simple and physically interesting Hamiltonians whose ground states are given by MPS. 
Clearly the symmetry breaking MPS is not the unique ground state of (lllj) . however 
the other ground states, can be found by applying the symmetry operators of su{2) 
group to the MPS. 

The other model that we find, labeled as model II, turns out to correspond to 
a family of solvable three-state vertex models on square lattices, [261 127] . In this 
model the degeneracy of the ground states, shows itself in a completely different way, 
namely, we find that the matrix product depends on a continuous parameter, but the 
parent Hamiltonian does not, i.e. 

H\^(g))=E \^{g)). (12) 

Thus if we expand the matrix product state \ip(g)) in terms of the parameter g, in 
the form 

N 

hA(<?)> = £<Wn>, (13) 
ra=0 

we obtain a large number of states \ip n ), which all have the same energy and thus 
represent part of the degenerate ground states of the Hamiltonian. In this way the 
MPS plays the role of a generating state for a set of degenerate ground states of the 
Hamiltonian, none of which has a MPS representation. The degree of degeneracy of 
the ground states increases with system size, and each state \t/j n ) has a complicated 
structure, and can not be represented as a matrix product, and thus the calculation of 
any of its correlation functions, or even its normalization, is quite difficult. However 
from the fact that the generating state \ip(g)) is a matrix product state, we can 
determine such correlations in closed form. 
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The structure of this paper is as follows: in section ([2]) we briefly review the matrix 
product formalism, with emphasis on the symmetry properties of the state and the 
parent Hamiltonian, in section ([3]) we consider three dimensional auxiliary matrices 
and classify them according to symmetries of the states which are constructed from 
them, namely symmetry with respect to rotation around the z-axis and the discrete 
parity and spin- flip symmetries. In this way we arrive at two specific forms of auxiliary 
matrices and consequently two specific models. Sections dH) and ([5]) are devoted to 
the detailed study of the above two models. The paper ends with a conclusion and 
an appendix. 



2 Symmetries of matrix product state and its parent 
Hamiltonian 

From (TjQ) we see that the collections of matrices {Ai} and {\UAiU~ 1 }, where A is a 
scalar, both define the same matrix product state. This freedom allows us to study 
the symmetries of MPS. A MPS will be symmetric under parity provided that we can 
find a matrix II such that 

n^ir 1 oc A T m 

and invariant under spin flip transformation, if we can find a matrix such that 

QA m Vl oc Ajn, 

where here and hereafter, A^z stands for j4_ m . As for continuous symmetries, consider 
a local symmetry operator R acting on a site as R\i) = Rji\j) where summation 
convention is being used. R is a d dimensional unitary representation of the symmetry. 
A global symmetry operator 1Z := R® N will then change this state to another matrix 
product state 

*W-i w — *' == HA' h A! i2 ■ ■ ■ A[ N ), (14) 

where 

A\ : RijAj. (15) 

A sufficient but not necessary condition for the state |\&) to be invariant under this 
symmetry is that there exist an operator U (R) such that 

R^Aj = U{R)A i U- 1 {R). (16) 

Thus R and U(R) are two unitary representations of the symmetry, respectively 
of dimensions d and D. In case that R is a continuous symmetry operator with 
generators T a , equation (|16p. leads to 

(T a ) ij A j = [T a ,A i ], (17) 

where T a and T a are the d— and D— dimensional representations of the Lie algebra of 
the symmetry. 



4 



A symmetric MPS need not be the ground state of a symmetric family of Hamilto- 
nians. To find the symmetric family of Hamiltonians we should construct the Hamil- 
tonian in a specific way. Let us first review how the Hamiltonian is constructed in 
general. From a matrix product state, the reduced density matrix of k consecutive 
sites is given by 



tr{{Al-..Al®A n . . ■ Aj k )E N ~ ! ' - 
' n '"' 1 tr(E N ) 



The null space of this reduced density matrix includes the solutions of the following 
system of equations 

d 

E c h-j k A jl ---A jk =0. (19) 

Given that the matrices A4 are of size DxD, there are D 2 equations with d k unknowns. 
Since there can be at most D 2 independent equations, there are at least d k — D 2 
solutions for this system of equations. Thus for the density matrix of k sites to have 
a null space it is sufficient that the following inequality holds 

d k > D 2 . (20) 

Let the null space of the reduced density matrix of k consecutive sites, denoted by 
Vfc, be spanned by the orthogonal vectors \e a ), (a = 1, • • • , s > d k — D 2 ). Then we 
can construct the local Hamiltonian acting on k consecutive sites as 

s 

h:=Y^ Ja\e a ){e a \, (21) 
a=l 

where J Q 's are positive constants. These parameters together with the parameters of 
the vectors \e a ) inherited from those of the original matrices Aj, determine the total 
number of coupling constants of the Hamiltonian. If we call the embedding of this 
local Hamiltonian into the sites I to I + k by /in+fc then the full Hamiltonian on the 
chain is written as 

N 

H = J2hl+k- (22) 
1=1 

The state is then a ground state of this Hamiltonian with vanishing energy. The 
reason is as follows: 

N 

(*\H\V) = tr(ff|*)(W|) = X>(>H.J+fcPi.J+fe) = 0. (23) 

l=i 

where pi t k+i is the reduced density matrix of sites I to I + k and in the last equality 
we have used the fact that h is constructed from the null eigenvectors of p for k 
consecutive sites. Given that H is a positive operator, this proves the assertion. For 
the Hamiltonian to have the symmetries of the ground state, the basis vectors of the 
null space should be chosen so that they transform into each other under the action of 
symmetries and the couplings J a should be chosen appropriately, see [28] for a more 
detailed discussion of this point. 
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3 Three dimensional auxiliary matrices 



As is clear from (|20p a sufficient condition for the existence of a null space V2 for a 
spin-one system is that the dimension of the matrices satisfy D 2 < 9 which restricts 
D to 1 and 2. The case of D = 1 has already been considered in [29J outside the 
framework of MPS formalism, and the case of D = 2 has been worked out in [2] and 
[5] as mentioned in the introduction. However we should emphasize that this is a 
sufficient and not a necessary condition and indeed we can take D > 3 and still find 
a non-empty null space V2, since the system of equations may not all be independent 
of each other. 



In this article we want to study in detail the case D = 3 and consider all the 
possible models which allow certain plausible symmetries, i.e. rotational symmetry 
around the z axis in spin space, and symmetry under parity and spin flip operations, 
these are the symmetries which have been taken into account in building optimal 
ground states for various models [31 13 El El El HH [HI [HI [22] . 

So let us take 3-dimensional matrices Ai, An and Aj and demand rotational sym- 
metry around the z axis in spin space. According to (|17p . this is equivalent to the 
following equations 



where S z 



[S t ,A 1 ]=A 1 , [S z ,A o ] = 0, [S Z ,A T ] 

diagonal(l, 0, —1). The immediate solution of these equations is 

/ 
c 
v d 




An 




(24) 



(25) 



where a,b,c,d,g,h and i are real parameters. By a transformation Aj — ► SAiS~ l 
where S = diagonal (1, a, ab) we can set the parameters of A\ equal to 1. Symmetry 
under parity now requires that there is a matrix II such that 



iL^rr 1 

A straightforward calculation gives 

n = 



A v 



and 



Ai 










1 





1 





1 
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h 
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t 











A T = 




c 










V 





c 






Finally we come to the symmetry under spin flip \m) 



\m) 



(26) 



(27) 



(28) 



It is readily seen 



that with these matrices the spin-flip symmetry is automatic, namely we have 



(29) 
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in which 

/ l\ 

n =\ o c o . (30) 

\c 2 J 

In order to guarantee that the matrix product state constructed in this way is 
the ground state of a Hamiltonian with nearest neighbor interaction, we consider the 
equation 

(HjAiA^O, (31) 

i,j=l,0,T 

which can be re-written as a matrix equation for the coefficients Cy in the form 

J2M kltljClj = 0. (32) 

i,j 

To have a solution we set 

det(M) = 0. (33) 
The determinant of M is readily calculated from its explicit form and is given by 

det(M) = (g 2 -h 2 ) 2 {2g 2 -h 2 )c 4 . (34) 

The vanishing of the determinant puts constrains on the parameters, namely we 
should have either c = 0, or h = ±g or h = ±y/2g, each choice leading to a different 
exactly solvable model. We omit the case c = since it leads to the condition A- = 
and hence reduces the model to an effectively two-state model, moreover in this case, 
spin- flip symmetry is lost due to the non-existence of an invertible matrix Q. Also 
it turns out that the other models with minus signs are equivalent to models with 
plus signs, see appendix A for a demonstration of this fact. So we are left with two 
different models which we label accordingly as model I (when h = V2g) and model 
II (when h = g) and study them separately in subsequent sections. 

Before proceeding to the models, we need to clarify a point about the number of 
parameters. Throughout our analysis we take N, the size of the lattice to be an even 
number. It appears that we have two continuous parameters in the matrix product 
states, namely g and c. None of these parameters can be gauged away by similarity 
transformations or scaling of the auxiliary matrices. However the MPS depends on 
only one parameter. To see this, let us expand the MPS in terms of the states \tp n ) 
which are defined to be linear superposition of all states which have n 0's. Note that 
for an even N, n will also have to be even. Since in the space of one site, the operators 
Ai and Aj act as raising and lowering operators, the trace of any string of operators 
^o^i an d A- is non- vanishing only if this string contains an equal number of A\ 

N — n 

and Aj. Thus any state \ip n ) comes with a coefficient g n c~^~ . Consequently for the 
un-normalized MPS we have 

N N 2 

mg,c)) = Y.9 n c^Wn) =cT £(-) 5 h/fe>- (35) 

n=0 n=0 C 



7 



Thus the normalized state and all the correlation functions will depend on only one 

2 

single parameter, namely For this reason we can put c = 1 and so the MPS will 
depend on only one single parameter g. 



4 Model I 

In this section we study in detail model I. The auxiliary matrices are 





( 


1 






( 


1 









(o 







Ai = 








S) 


, A = g 









i) 









i) 




v° 









\ 














1 





(36) 



First we derive the one-parameter family of parent Hamiltonians and then calculate 
the one and the two-point functions. 

4.1 The Hamiltonian 

Here we have h = V2g and the null space V2 is spanned by one single vector 

|e) = |00)- 5 2 |l,T)- 5 2 |T,l). (37) 

We take the local Hamiltonian as hj = n^^p l e )( e l from which the full Hamiltonian 
turns out to be 

ffj(u) = (l-u 2 )N+J2 Sl l Sl i+1 -2(l+2u)Sl i +u 2 (S l -S i+1 ) 2 +u{S i -S l+1 , S z>i S z>i+1 }, 

(38) 

2 

where u := jz^z- When g = u = 0, the Hamiltonian turns into 

^i = E(^- 1 )(^i+i- 1 )- ( 39 ) 

i 

In this limit, since Aq — > 0, the MPS becomes an expansion of states consisting only of 
1 and 1. Such a state is clearly the ground state of the Hamiltonian (|39p . however, the 
Hamiltonian ([39]) has a highly degenerate ground state, which is not captured by the 
MPS in this limit. In fact, any basis state in which no two 0's are adjacent is a ground 
state of this Hamiltonian, with energy Eq = 0. The number of ground states of H\ is 
equal to tr(A ), where iV is the system size and A is the adjacency matrix in which 
allowed adjacent configurations are designated by 1 and disallowed configurations by 
0. In Hi the only configuration which lifts the local energy from to 1 is that of two 

/ill 

adjacent zeros, so in the basis which we have chosen, A = 



V 



1 1 I . For large 
1 1 1 



N, we will have tr(A) N ~ (1 + y/3) N , where 1 + y/3 denote the largest eigenvalue of A. 
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When g — > 1 (or u — > oo), Hj(u) turns, modulo a multiplicative coefficient, into 
the following Hamiltonian, 

ff 2 = X)[(S < -S i+ i) 2 -l]. (40) 

i 

This is a simple and interesting Hamiltonian and thus our result implies that its 
ground state is of the form of an MPS, with matrices given by (|36h . for g = 1. 

Note that in the limit g = 1, the null eigenvector (|37|) becomes a singlet, the 
spin-0 representation of angular momentum, implying that the Hamiltonian should 
be a scalar which conforms with the form of the Hamiltonian (|40|) . However the 
matrices (|36p (for g = 1) do not transform as a spherical tensor operator under 
angular momentum, that is the following relations are not satisfied as required from 
(HZ]): 

[Sj^-Ajjj] = mA m 

[S+,A m ] = J2-m(m+l)A m+1 



[S-,A m ] = v /2-m(m-l)A m _ 1 , (41) 

where S z , and S± are the three dimensional representation of su{2). This means that 
the Hamiltonian is symmetric under the full rotation group, but the ground state, 
breaks this symmetry. Therefore other degenerate ground states can be constructed 
by applying rotation operators to this state. However the actual degeneracy is much 
larger and it grows exponentially with system size as 

(*±fi)N j3o]. 



We should stress that the Hamiltonian H% corresponds to a particular point in the 
class of bilinear-biquadratic spin-1 chains with the Hamiltonian 

N 

H{6) := £ [cosO Si ■ S i+1 + S m9 (Si ■ S m ) 2 ] (42) 

8=1 

which has extensively been studied by various methods |314 13"2"l I3"3"l I3"4"l [331 136] . 
4.2 Correlation functions 

The correlation functions of this model are determined after lengthy but straightfor- 
ward calculations starting from © and ([6])- In the thermodynamic limit, the results 
are as follows: 

(S z ,i) = (S^) = 0. (43) 

Thus in the ground state, there is no magnetization. Note that due to rotational 
invariance in the x — y plane of spin space, in all the correlation functions below, we 
can change x to y or any other direction in the x — y plane. To describe the other 
correlation functions, let us introduce the parameter 7 := ^/ g 4 + 8. Then we have 

(S* ) = 8 - (S 2 •) = g2(37 + ff2) + 4 (44) 

{z ' l) 7(35 2 + 7)' {x ' %) 7(35 2 +7) ' 1 j 
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g 

Figure 1: The average values (Sf) and (S^) for model I. In the limit g — > oo(u — > — 1) 

, the spins lie in the x — y plane. 



Figure ([T]) shows the plot of (<Sj ^) and (<Sji) as a function of the parameter g. In 
the limit g — > 0, we have (S^i) — * 1) implying that there is no in the expansion of 
the state. This is in accord with our picture of the MPS, since in this limit Aq — > 0. 
Finally for the two point correlations of longitudinal and transverse components of 
spins we find 

(S x ,xS x , r ) = G ± (g)e-^/^, (45) 
where the magnitude of correlations are given by 

4(g 2 + 7 ) r , v 8ff 2 (g 2 + ^ + 7 ) 2 (Af , 



and the longitudinal and transverse correlations are given by 

1 

_ & 



These are plotted in figure ((2J). Note that G||(y) = (S Z) iS z ,2) an d G±(g) = {S x> iS x ,2}- 
In the limit 5 — > 0, where Hj becomes an Ising-like Hamiltonian, the above equations 
show that transverse correlations vanish, and longitudinal correlations approach the 
value —1/2. However we should note that in this limit, the ground state is highly 
degenerate. In fact as stated above, any basis state in which there are no two adjacent 
0's is a ground state of H\. However the MPS does not capture this degeneracy, but 
is only one of the many ground states. In the limit g — * 1, where the Hamiltonian 
turns into (j4T)]) . we have (5* 2 j) = | and (5* 2 j) = I, and the correlation lengths tend 

t° *|| = R3J = - 910 ^ = In(3/(1W2)) = 4603 - 
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-0.2- 



-0.4. 



a b 

Figure 2: The correlation functions for longitudinal and transverse components of 
spins, in model I; a) Correlation lengths, b) the magnitudes of correlations. 



5 Model II 

For this model, the auxiliary matrices are 

An-- 





5.1 The Hamiltonian 

The null space V 2 is spanned by the following two vectors: 

1 




(48) 



|ei; 

\ e 2) 



V2 
1 

72 



(]01> - |10», 

(|ol) - |T0». 



(49) 



These vectors are eigenvectors of the local two-site S z operator, are invariant under 
parity and transform into each other under spin flip transformation. Therefore if we 
take the local symmetric Hamiltonian as 



hn = |ei)(ei| + |e 2 )(e 2 |, 



(50) 



where we have set a total multiplicative constant equal to unity, the final total Hamil- 
tonian is spin-flip and parity invariant and moreover commutes with the third com- 
ponent of spin, i.e. [hji, S z ] = 0. Its explicit form in terms of local spin operators can 
be determined after some algebra: 



N 



Hn — ^ ^Sgj — {Sj.Sj+i, S Zj iS z> i+\} — Sj.Sj+i + S Zj iS Zj i+i. 



(51) 



»=i 
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1 

a 





b 



|a-cl= b 



1 

c 



o 

a 



Figure 3: The independent Boltzmann weights of the two dimensinal vertex model 
corresponding to the quantum spin chain (|5ip . 



This Hamiltonian was first found in [27] . The history is the following: The exhaustive 
solutions of the Yang-Baxter equation, corresponding to a three-state 19- vertex model 
on a square lattice were first found in [26]. These solutions reproduced many of the 
already known exactly solvable vertex models in addition to four new models. These 
models, labeled I, II, III and IV in [27], were then studied in detail in [27] ■ where 
the thermodynamic properties of these new models, including the partition function 
and correlation lengths were derived. Two of these models, namely models I and II, 
however allowed a more complete solution (due to the so called crossing symmetry 
of the R matrix, the solution of the Yang-Baxter equation) which allowed the exact 
determination of the ground state energy per site. However the other two models, 
models III and IV, lacked this symmetry, and no exact solution for the ground state 
energy was given. It could however be established that such models can be mapped to 
6- vertex models, i.e. two state models with 6 allowed configurations. The Hamiltonian 
(|51h corresponds in fact to the Hamiltonian of model III in [27], for A = 1, where A is 
a particular combination of Boltzmann weights. For its definition and the Boltzmann 
weights see j27J. According to [27], the three-state 19- vertex model corresponding to 
this spin chain is defined by the Boltzmann weights shown in figure ([3]). 

Note that each Boltzmann weight R^'g corresponds to the spin labels, /i, v, a, (3 
on left, right, bottom and top links of a vertex and the R matrix and correspondingly 
the Boltzmann weights have the following symmetries (with a := —a, etc) : 

= RZ = R V A = fl?J. (52) 

The exact correspondence between the Hamiltonian ()51|) and that of model III of 
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|27j is as follows: 



Hn — ^6-vcrtex + N, 



(53) 



where N is the number of sites of the chain. This correspondence allows to determine 
exactly the ground state energy of the corresponding 6-vertex model, for the range 
of values of Boltzmann weights of the 6-vertex model shown in figure ([3]). Since the 
ground state energy of H is zero, we find the ground state energy of ii~6-vertex to be 
Eq = —N, giving an energy per site equal to eo := = — 1. 



5.2 The explicit form of the ground states 

From the form of the local Hamiltonian one finds that any state which comprises of 
only l's and l's, with no O's, is a ground state of this Hamiltonian. The number 
of such states is 2 N , where N is the size of the system and their explicit form is 
|si, S2, • • • s^), where s, £ {1, 1}. However there are other non-trivial ground states, 
those which contain also a number of O's. The interesting point is that these kinds of 
ground states, are nicely captured by the MPS ((151) . We can find ./V such degenerate 
states. To see this we note that MPS ([35]) depends on a continuous parameter g while 
the Hamiltonian does not. This dependence is genuine and can not be gauged away 
by similarity transformation A m — ► U A m U~ l . Let us expand the MPS in terms of 
powers of g in the form 

N 

mg)) = Y,9 n \*Pn), (54) 

n=0 

where 

IV>n) = tr(A h A i2 ■■■A iN )\ii,i 2 ,---iN), (55) 

h,i2,---iN 

and implies that in each term only n Aq^s exist. 

In view of (|48j) and the definition of MPS (JH) the powers of g enumerate the num- 
ber of O's in each state and so each \ip n ) is a superposition of states each of which has 
exactly n local O's and an equal number of l's and l's. Note that since we have taken 
./V to be even, this implies that n is also an even number. 

To determine the explicit form of each state \tp n ), consider the form of the matrices 
Aq, A x and Aj in ([18]). With 




we have 

Ai = |l)(0| + |0)(T|, A = gl, A-=|T)(0| + |0)(1|, (56) 
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N-4 

i<j<k<l 



Figure 4: The state |i/^iv— 4) • All the dots are at the state |0), while those which are 
connected are at the state |1,1) + |1,1). 

or more compactly 

A m = \m)(0\ + \0)(m\, A =gl, (57) 

where m = 1, 1 and / is the identity matrix. The product of any string of matrices 
A m , m = 1, 1 has a simple structure. One can easily show that 

3rri2 ,m3 ^4 ,1715 ' ' ' &rri2K-2,m2K-i l^l) (^2fc | ■ V*°) 

This allows us to determine the explicit form of any of the states \ip n )- Consider 
for example \iPn), where there is no 1 or 1. This is a simple state \iPn) = 3|0, 0, • • • 0), 
where the factor 3 comes from taking the trace of the identity matrix, related to the 
product of all the Aq matrices. The next state is \1pN-2} in which two of the zeros 
have been replaces by 1 and 1. Using (j58H we find that this is a state of the form 



\if> N - 2 ) =2 ^ (|0, • • • 0, l m , 0, • • • , 0, 0, • • • 0> + |0, • • • 0, l m , 0, • • • , 0, 1„, 0, • • • 0)), 

(59) 

where the two non-zero spins occur at sites m and n respectively. Let us define 
\aij) := lj) + lj) where the indices denote the sites of the lattice, and all other 
sites are occupied by 0. Then I^A^— 2) = %J2i<j \ a ij)- 111 view of ([58]) . we find the 
structure of all other states \if) n ). For example we have 

\ipN-i) = (\ a ij)\ a kl) + \otjkOtli), (60) 

i<j<k<l 

which is pictorially depicted in figure ([!]) and 

\1pN-6) = (\^)ij\a) k i\a) mn + \a) n i\a) jk \a)i m ) . (61) 

i<j <k<l<m<n 

Finally we find v/jq, in which there are no zeros, and has a dimmerized or Majumdar- 
Ghosh like structure, namely 

\ipo) = I«i2)|a34) • • • \&N-i,n) + |"23)|a45) • • • |«iV,i), (62) 
which is depicted in figure (0). 
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¥ = 





Figure 5: The state \i/jq) has a Majumdar- Ghosh like structure. Note that a line here 
means only an entangled state |1, 1) + |1, 1) and not a spin singlet. 

5.3 Correlation functions 

We have found that the number of degenerate ground states of model II, for a chain of 
size N, is at least 2 N + N. Of these, the 2^ states \s\, S2, • • • sn), where Sj G {1, 1} are 
uncorrelated, even for finite-size systems. The other N states are non-trivial and have 
MP representations. Here we calculate the correlation functions for the other types 
of states which are correlated. Were it not for the matrix product formalism, such 
calculation could be very difficult. As a first step let us determine the normalization 
of the states \ip n ), and then proceed to the calculation of correlation functions of 
various operators. 

Let us first fix our notation for matrix product operators, using their general 
definition (jSJ). We have 

E Sz = Ai ® A x - A T <g> A- =: U, 
E S 2 = Ai (g> A ± + A T <g> A T = V, 
E = Ai ® A x + A T (8 A- + A A = V + g 2 I. (63) 

Let us start by finding the normalization of the states \ip n )- It is obvious that for 
m ^ n, (ipmltpn) = ; since these two states, have different number of zeros in their 
expansion. Therefore 

N 

WgM(g)) = Y,9 2n (^n). (64) 

n=0 

The left hand side is found from ([2]) to be 



MgM(g)) = He n ) = Hv + g 2 i) N = T J g 2n [ J tr(v N ~ n ). 




Comparing with the previous formula, we find 

(iMVfe) = ( n ) HV N ~ n ). (65) 
5.3.1 One-point functions 

In the thermodynamic limit iV — ► 00 (n=finite), the traces are simplified by taking 
only the largest eigenvalues, but before going to this limit, let us derive closed formulas 
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for some of the correlation functions. Prom the nature of the states \ip n ), it is obvious 
that 

{lp n \S z \lp n ) = {lpn\S a \lp n ) = 0, (66) 

where a is any direction perpendicular to the z axis. To obtain the average of S 2 , we 
note again that for in ^ n, (Tp n \S 2 \ip m ) = and hence 

N 

^{g)\Sl\i>(g)) = £ g 2n (^ n \S 2 \^ n ). (67) 

n=0 

The left hand side is obtained from ©: 

N ~ l ( N 1 \ 



(68) 



n=0 

which in view of (|67|) gives 



N-l\ 

(^ n \S 2 ^ n ) = V " ) = N-n 



(69) 



as expected, since the operator counts the average number of l's or l's in the state. 
From (S 2 + S 2 + Sf) = 2 and the rotational symmetry around the z axis, we find for 
any unit vector a in the x — y plane 

&n\S 2 M = N + n 

5.3.2 Two-point functions 

Let us now derive the two point correlation functions, {tp n \S Zt iS Zt r\tp n ) . We have 

(^(g)\Sz,iS z , r \i>(g)) = tr(E Sz E r - 2 E Sz E N - r ) = tr(U(V + g 2 I) r - 2 U(V + g 2 I) N - r ) 

Expanding the right hand side and using the diagonal property of this two-point 
function, namely that for m ^ n, {ipn\S z ,iS z ^ r \ipm) = 0, we find 



{i>n\S z ,lS Z) r\i>r, 



E *=° ( r k 2 J ( n - k ) HUVr- 2 -*UV^-n +k) 



(ipnl^n) ( N 



(71) 



n 



tr (yN- 



The expansion of {ij){g)\S 2 1 S 2 r .\ijj(g)) gives a simple result, since in this case, the 
matrix product operator E S 2 = V commutes with E, that is 
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(V^I^XrlV'G?)) = tr{V(V + g 2 ir 2 V(V + g 2 I) N - r ) = tr{V\V + g 2 I) N ~ 2 ) 
= £ ( N ~ 2 ) 9 2n HV»-% (72) 

giving the final distance-independent result 

( N ~ 2 ) 

mS 2 Ztl Sl r \^ n ) = \ n J = (N- n)(N - n - 1) 
(V-nlVn) / N \ N(N-l) 



(73) 



Finally we come to the two-point function (S x i8 xr ) between transverse compo- 
nents of spins. Here we encounter an essential difference with the previous cases, in 
that the operator S x is not diagonal between the states \ip n ), i.e. {il) n \Sx,iSx,r\' t l>m} 
for m 0. To circumvent this problem, we do the following expansion, for g and h 
two arbitrary variables: 

TV 

(ip{g)\Sx,iS x ,r\ilf(h)) = 9 n h m mS xA S x , r \ij m ). (74) 

n,m=0 

To calculate the left hand side within the matrix product formalism, we should slightly 
adapt equations ((2J), © and © to conform to the present situation. Let \4>) and \ tp) be 
two matrix product states on the same periodic lattice of size N, defined by matrices 
{Ai} and {Bi} respectively. Then by following the steps which led to ([2]) and ([5]) we 
find 

<<f#> =tr(E N ), (75) 
where E = J2i A* <8> Bi. The one point function of an observable O will be given by 

<0|O|V> = tr{E E N - x ), (76) 

where 

E = Y f A$®B j (i\0\j). (77) 

»>j 

We now use the above formulas to calculate the left hand side of (|74p . where \ip(g)) 
and \ip(h)) are both matrix product states with matrices Ai(g) and Ai(h) respectively. 
For the operator in question, namely S x , we find from 

5 X = -^(|0)(1| + |0)(T| + |1)(0| + |I)(0|) 

the form of its matrix operator, 

E Sx = ^=(gI®(A l +A-) + {A l + A-) ® hi) , (78) 
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which can be simply written as 

E Sx = ^=(gX 2 + hX 1 ), (79) 

where X = A\ + Aj and the indices on X means its embedding on the first and second 
spaces. Thus we obtain 

(iP{g)\S Xil S x , r \Tj;{h)) = \tr{(gX 2 + hX x )(V + ghI) r ~ 2 (gX 2 + hX 1 )(V + gh) N ~ r ~ 2 ). 

(80) 

Expanding the right hand side and collecting the coefficient of {gh) n will give 



(lpn\S X) lS Xyr \^ n ) _ 1 



£C5ri 2 H:_-r_-* 2 W 



L N,r,k,n> 



(^n\tpn) 2 I N 

\ n 



(81) 



tr(V N ~ n ) 



where 



^N,r,k,n := X 2 V r ~ 2 - k X 1 V N - r '- n - 1+k + X 1 F r - 2 - fc X 2 F iV ~ r - n ~ 1+fe (82) 
We now come to the thermodynamic limit of these correlation functions. 

5.3.3 The thermodynamic limit 

In this subsection we calculate the thermodynamic limit of the correlation functions 
found above, by which we mean taking N — > oo while keeping r and n fixed. In this 
limit, the sums in the numerators of equations (|7ip or (|8ip are dominated by their 
k = term, therefore we find from (|71|) 

(^nlS^S^n) = tr(UV r ~ 2 UV N - r ~ n ) 

(^n|Vn) tr{V N ~ n ) ' 1 ' 

Moreover when taking the traces of V N , we need only keep the eigenvalue with largest 
magnitude. The matrix V has two eigenvalues with largest magnitudes, namely A± = 
±\/2, whose corresponding eigenvectors we denote by |±). Their explicit form is 

|±) = i(|ll)±v^|00> + |Tl)) > \± = ±V2. (84) 

Therefore by using the fact that N and n are even (see the discussion preceding 
equation (|35j) ). we arrive at 

{il>n\s z ,is z , r \ii> n ) _ i {+\uv r - 2 u\+) + {-iy {-\uv r ' 2 u\-) 

(VnlV-n) 2 2-/2 " ^ > 

Using equations (j63[) and (|84l) and noting the explicit form of the matrices Aq = 
|1)(0| + |0)(1| and A\ = |0)(1| + |1)(0|, we find after a straightforward calculation, that 
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the above matrix elements vanish unless r = 2. Using the fact that (±|C/ 2 |±) = — 1, 
this will give the correlation function 

, ; r = -qdr.2- (86) 

Similarly for the correlation of transverse components we find from (|8ip and following 
a similar reasoning as above that 



(87) 



(i> n \i> n ) 2Ntr{V*-n)i 

which simplifies to 

{lpn\Vn) N 2 — 

where VL X = X\V r ~ 2 X2 + X^V" X\- Obviously this tends to zero in the limit we 
are considering, due to the pre- factor jjj. 



6 Discussion 



We have made a detailed study of two new spin systems whose ground states can be 
found exactly within the matrix product formalism, and have shown that the method 
of MPS formalism may be more fruitful than it appears at first sight. In the space of 
quantum spin chains solvable by this formalism, there may be physically interesting 
models with rich properties. For example, there may be ground states which break the 
continuous symmetries of their parent Hamiltonian, or one may find matrix product 
states which capture as a generating state, a large number of degenerate ground 
states of a given Hamiltonian. In this paper we have made an exhaustive study 
of spin-1 matrix product states with three-dimensional auxiliary matrices and have 
found two distinct class of models. The first model is a one-parameter family of 
models interpolating between the Ising-like Hamiltonnian 

i 

and the rotationally invariant Hamiltonian 

H2 = y^(Sj • Sj + i) 2 . 

i 

The ground state of this latter Hamiltonian which is in the form of an MPS, breaks 
the rotational symmetry SO(3) of H2 to that of SO(2). Symmetry generators now will 
give other degenerate ground states of H2. The second model which we have found 
corresponds to a kind of 6- vertex model, already studied in [27]. Here the MPS is a 
generating state of the degenerate ground states of the Hamiltonian, none of which 
may have a simple MPS representation. By suitable manipulations of this generating 
state, we have been able to find one and two-point functions for these ground states, 
which are other-wise very difficult to calculate. 
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8 Appendix A 

In this appendix we show explicitly for model II, that the two different choices h = g 
and h = —g lead to equivalent models. A similar analysis applies to model I, which 
we omit here. Let h = ag, where a = ±1. Then it can be verified that the null space 
V2 is spanned by the vectors: 

| ei ) = i(|01)-a|10)), 

|e 2 > = ^(|01>- < 7|T0)). (89) 

Writing the local Hamiltonian as usual tin = |ei)(ei| + |e2) (e2 1, we see that 

huip = -1) = (I ® R z (7r))h n (a = 1)(J ® ^(vr)- 1 ), (90) 
where R z (ir) is a rotation around z axis by tt degrees. Thus we arrive at 

H II (-a) = KH n (a)K- 1 , (91) 
where we have taken N to be an even number and 1Z is a global rotation of the type 

K = (g)R? +1 (ir). (92) 

i=0 

This shows that Hu(a) and Hjj(-a) are iso-spectral and have the same thermody- 
namic properties and for this reason we have considered only the model Hu{a = 1). 
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